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ABSTRACT 

The satellite-borne experiment PAMELA has been used to make new measurements of cosmic ray H and He isotopes. 
The isotopic composition -was measured bet-ween 100 and 600 MeV/n for hydrogen and bet^ween 100 and 900 McV/n 
for helium isotopes over the 23'^'^ solar minimum from July 2006 to December 2007. The energy spectrum of these com- 
ponents carries fundamental information regarding the propagation of cosmic rays in the galaxy -which are competitive 
with those obtained from other secondary to primary measurements such as B/C. 

Keywords: Astroparticlc physics, cosmic rays 



1. INTRODUCTION 

Hydrogen and helium isotopes in cosmic rays are gen- 
erally believed to be of secondary origin, resulting from 
the nuclear interactions of primary cosmic-ray protons 
and *He with the interstellar medium, mainly through 
spallation of primary ^He nuclei or through the reac- 



tion p + p 



'U 



7r+. These isotopes can be used 



to study and constrain parameters i n propagation mod- 
els for galactic cosmic rays (GCRs) (jStrong et al. 1 120071 : 
lTomassettill2012t ICoste et aLll2012[ ). ^h and ^He are 
the most abundant secondary isotopes in galactic cosmic- 
rays and have peculiar features: ^H is the only sec- 
ondary species (apart from antiprotons) that can also be 
produced in proton-proton interactions and ^He is the 
only secondary fragment with a A/Z significantly differ- 
ent from two {A and Z being respectively the mass and 
charge number). 



The importance of light isotopes has been known for 
about 40 years, when the first meas u rements became 
available (IGarcia-Munoz et al. |[T975ai rbl: [Mewaldt et al. I 
119761: ILeech fc O' Gallagheri I1978D . Measurements re- 
quire very good mass resolution, a challenge for in- 
struments deployed in spac e. With the excep tion of 
the results from AMS-01 (jAguilar et al. I |2011[ ) most 
of the measure ments were perfo rmed using strato 
spheric balloons (iWang et al. II2002I: JReimer et al. 199"! 
Wefel et al. I 119951: iWebber et al. I Il991i: iBeattv et al. 



1993D . where the residual atmosphere above the instru- 



ment caused a non-negligible background of secondary 
particles. The atmospheric background estimation is 
subject to large uncertainties (e.g. the limited knowl- 
edge of isotope production cross sections). Due to 
these limitations, experimental errors are generally very 
large and the focus of measurements therefore shifted 



to other secondary species, l ike boron or sub-iron nuclei 
(|Strong fc Moskalenkol[T998l) . 

The light-isotope quartet offers an independent unique 
way to address the issue of "universality in GCR propa- 
gation" , which can be crucial when analysing antiproton 
and positro n spectra to sear ch for possible primary sig- 
nals (see e.g. lStepheni (|1989[ )V The PAMELA experiment 
has been observing GCRs over the 23'^'^ solar miminum 
since July 2006 at an altitude ranging from ~ 350 km to 
^ 600 km on-board of the Russian Resurs-DKl satellite 
which executes a quasi-polar orbit. The low-earth orbit 
allows PAMELA to perform measurements in an environ- 
ment free from the background induced by atmospheric 
cosmic rays. 

The results presented here are based on the data set 
collected by PAMELA between July 2006 and Decem- 
ber 2007. From about 10^ triggered events, accumulated 
during a total acquisition time of 528 days, 5,378,795 
hydrogen nuclei were selected in the energy interval be- 
tween 100 and 600 MeV/n and 1,749,964 helium nuclei 
between 100 and 900 MeV/n. The data presented here 
re place and complete th e preliminary results presented 
in iCasolino et alA (|2011l ). A more complete evaluation 
of the selection efficiencies and contamination due to an 
improved simulation resulted in better reconstructed ^H 
and '^He spectra and '^He/^He ratio. 

2. THE PAMELA APPARATUS 

The PAMELA spectrometer was designed and built to 
study the antimatter component of cosmic rays from tens 
of MeV up to hundreds of GeV and with a significant in- 
crease in statistics with respect to previous experiments. 
To reach this goal the apparatus was optimized for the 
study oi Z = I particles and to reach a high level of 
electron-proton discrimination. 

The core of the instrument (Fig. [T]) is a permanent 
magnet with an almost uniform magnetic field inside the 
magnetic cavity which houses six planes of double-sided 
silicon microstrip detectors to measure the trajectory of 
incoming particles. The spatial resolution is '--^ 3 /im in 
the bending view (also referred as x-view) and ^^ 11 fim 
in the non-bending view (also referred as j/-view). The 
main task of the magnetic spectrometer is to measure 
particle rigidity p = pc/Ze (p and Ze being respectively 
the particle momentum and charge, and c the speed of 
light) and ionization energy losses (dE/dx). 

The Timc-of-Flight (ToF) system comprises three dou- 
ble layers of plastic scintillator paddles (SI, S2, and S3, 
as shown in Fig [T|) with the first two placed above and 
the third immediately below the magnetic spectrometer. 
The ToF system provides 12 independent measurements 
of the particle velocity, /3 = v/c, combining the time of 
passage information with the track length derived from 
the magnetic spectrometer. By measuring the particle 
velocity the ToF system discriminates between down- 
going particles and up-going splash albedo particles thus 
enabling the spectrometer to establish the sign of the 
particle charge. The ToF system also provides 6 inde- 
pendent dE/dx measurements, one for each scintillator 
plane. 

A silicon-tungsten electromagnetic sampling calorime- 
ter made of 44 single-sided silicon microstrip detectors 
interleaved with 22 plates of tungsten absorber (for a to- 
tal of 16.3 Xq) mounted below the spectrometer is used 
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Figure 1. A schematic overview of tiie PAMELA satellite exper- 
iment. The experiment stands ~1.3 m high and, from top to bot- 
tom, consists of a time-of-fiight (ToF) system (SI, S2, S3 scintil- 
lator planes), an anticoincidence shield system, a permanent mag- 
net spectrometer (the magnetic field runs in the y-direction), a 
silicon-tungsten electromagnetic calorimeter, a shower tail scintil- 
lator (S4), and a neutron detector. 

for hadron/lepton separation with a shower tail catcher 
scintillator (S4). A neutron detector at the bottom of 
the apparatus helps to increase this separation. 

The anticoincidence (AC) system comprises 4 scintil- 
lators surrounding the magnet (CAS), one surrounding 
the cavity entrance (CAT) and 4 scintillators surround- 
ing the volume between SI and S2 (CARD). The system 
is used to reject events where the presence of secondary 
particles generates a false trigger or the primary particle 
suffers an inelastic interaction. 

The readout electronics, the interfaces with the CPU 
and all primary and secondary power supplies are housed 
around the detectors. The apparatus is enclosed in a 
pressurized container attached to the side of the Resurs- 
DKl satellite. The total weight of PAMELA is 470 kg 
while the power consumption is 355 W. A more detailed 
description of the instruments and the data handling can 
be found in lPicozza et al. I (|2007l ). 

3. DATA ANALYSIS 

3.1. Event selection 

The month of December 2006 was discarded to avoid 
possible biases from the solar particle events that took 
place during the 13*^^ and 14*^^ of December. The event 
selections adopted are similar to those used in previ- 
ous works on the hig h energy proton and helium fluxes 
(jAdriani et al. |[20Tll ) a nd on the time depeii dence of the 
low energy proton flux (jAdriani et al. II2013I ). 

3.1.1. Event quality selections 

In order to ensure a reliable event reconstruction a set 
of basic criteria was developed. The aim of these require- 



ments was to select positively charged particles with a 
precise measurement of the absolute value of the parti- 
cle rigidity and velocity. Furthermore, events with more 
than one track, likely to be products of hadronic inter- 
actions occurring in the top part of the apparatus, were 
rejected. 
Events were selected requiring: 

• A single track fitted within the spectrometer fidu- 
cial volume where the reconstructed track is at 
least 1.5 mm away from the magnet walls. 

• A positive value for the reconstructed track curva- 
ture. 

• Selected tracks must have at least 4 hits on the x- 
view and at least 3 hits on the y-view to ensure a 
good rigidity reconstruction. 

• A maximum of one hit paddle in the two top planes 
of the ToF system. 

• The hit paddles in SI and S2 must match the ex- 
trapolated trajectory from the spectrometer. 

• A positive value for the measured time of flight. 
This selection ensures that the particle enters 
PAMELA from above. 

• For the selection of the hydrogen sample no ac- 
tivity in the CARD and CAT scintillators of the 
anticoincidence system is required. 

The anticoincidence selections on the hydrogen sample 
were necessary since most secondary particles that en- 
tered the PAMELA fiducial acceptance were by-products 
of hadronic interactions taking place in the aluminium 
dome or in the SI and S2 scintillators. Such particles 
were generally accompanied by other secondary particles 
which hit the anticoincidence detectors. For the selec- 
tion of the helium sample there were no anticoincidence 
requirements since contamination by secondary helium 
coming from heavier nuclei spallation is negligible. 

3.1.2. Galactic particle selection 

The Resurs-DKl satellite orbital information was used 
to estimate the local ge omagnetic cutoff, G, in the 
Stormer approximation (iShea et al. I Il987f) using the 
IGRF magnetic field model (iMacMillan fc Maui l2005l) 
along the orbit. The maximum zenith angle for events 
entering the PAMELA acceptance was 24 degrees with a 
mean value of 10 degrees. To select the primary (galac- 
tic) cosmic ray component particles were binned by re- 
quiring that pm > k ■ G, where pm is the lowest edge 
of the rigidity interval and fc = 1.3 is a safety factor 
required to remove any directionality effects due to the 
Earth's penumbral regions. Galactic particles losing en- 
ergy while crossing the detector may be rejected by this 
selection. This effect is accounted for using Monte Carlo 
simulations. 
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Figure 2. Energy loss in the silicon detectors of the tracking sys- 
tem (mean energy deposition in all planes hit) as a function of 
reconstructed rigidity for positively charged particles. The hydro- 
gen and helium bands are clearly visible. The black lines represent 
the selection for H and He nuclei. 
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Figure 3. Mass separation for Z = I {top) and Z = 2 {bottom) 
particles using the /3-rigidity method. The black lines represent the 
expectations for each isotope. 



12 dE/dx measurements. The arithmetic mean of those 
measurements is shown in Fig. [2l A rigidity dependent 
selection on the mean dE/dx from the spectrometer is 
used to select Z = 1 or Z = 2 candidates and is depicted 
by the solid lines in Fig. [2l 

The residual contamination of Z = 2 particles in the 
Z = \ sample was studied selecting helium events using 
the dE/dx information from the ToF system, and then 
applying the Z = \ selection previously described. The 
fraction of misidentified helium events was found to be 
less than 10""*. Since the ^H/^'^He ratio is roughly 0.15 the 
resulting contamination in the ^H sample from misiden- 
tified "^He was estimated less than 10^"^. Similarly the 
Z = \ residual contamination as well as contamination 
by heavier nuclei in the Z = 2 sample was estimated to 
be negligible. 



3.1.3. Charge selection 

Particle charge identification relies on the ionization 
measurements provided by the magnetic spectrometer. 
Depending on the number of hit planes there can be up to 



3.2. Isotope separation 

The selection criteria described in the previous section 
provided clean samples of Z = 1 and Z = 2 particles. 
In each sample an isotopic separation at fixed rigidity is 



possible by reconstructing /3, where 



P 



Z^p^ 



-1/2 



(1) 



as shown in Fig. [3] for events in the (/3, p) plane. 

Isotope separation as well as the determination of iso- 
tope fluxes was performed in intervals of kinetic energy 
per nucleon. Since the magnetic spectrometer measures 
the rigidity of particles, this implies different rigidity in- 
tervals according to the isotope under study. For exam- 
ple Fig. |4] shows the 1//3 distributions used to select ^H 
(top panel) and ^H (bottom panel) in the kinetic energy 
interval 0.329 - 0.361 GeV/n corresponding to 0.85 - 0.9 
GV for ^H and 1.7 - 1.8 GV for ^H. Particle counts were 
subsequently extracted from a Gaussian fit to the 1//3 
distribution in each rigidity range as shown by the solid 
lines in Fig. [4l 

Separation between '^He and ^He was obtained in a 
similar way. Fig. [5] shows the 1//3 distributions used to 
select ^He (bottom panel) and '''He (top panel) in the ki- 
netic energy interval 0.312 - 0.350 GeV/n corresponding 
to 1.24 - 1.32 GV for ^He and 1.65 - 1.76 GV for ^He. 

It should be noted that, because of the large proton 
background, an additional selection, based on the lowest 
energy release among the 12 measurements provided by 
the tracking system (often referred as truncated mean), 
was used to produce the 1//3 distributions in the ^H case. 
Fig. |6] shows this quantity for Z = 1 particles. The solid 
line indicates the condition on the minimum energy re- 
lease used for the selection. 

The selected number of hydrogen and helium events are 
summarized in the second and third column of Tables [1] 
andU 

3.3. Flux determination 

The procedure described in the previous section was 
used to estimate the number of ^H and ^H events in the 
Z = 1 sample and the number of "^He and ^He events 
in the Z = 2 sample. To derive each isotope flux the 
number of selected events had to be corrected for the 
selections efficiencies, particle losses, contamination and 
energy losses. These corrections were obtained using a 
Monte Carlo simulati on of the PAMELA appar atus based 
on the GEANT4 code (|Agostinelh et al. I [20031 ) and from 
the flight data. The simulation contains an accurate 
representation of the geometry and performance of the 
PAMELA detectors. The measured noise of each silicon 
plane of the spectrometer and its performance variations 
over the duration of the measurement were accounted for. 
The simulation code was validated by comparing the dis- 
tributions of several significant variables with those ob- 
tained from real data. Hadronic interactions for all the 
isotopes under study were handled via the QGSP_BIC_HP 
physics list. 

The following corrections to the number of selected 
events were applied: 

• Selection efficiencies: The redundant information 
provided by PAMELA allowed most of the selec- 
tion efficiencies to be estimated directly from flight 
data. For example the efficiency of the charge selec- 
tions was evaluated on a sample of events selected 
with the ToF dE/dx measurements in the same 
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Figure 4. Examples of 1//3 distributions for iiydrogcn in the 0.329 
- 0.361 GeV/n kinetic energy range for ^H (top) and ^H (bottom). 
The solid line shows the estimated ^H and ^H signal while the 
dashed line shows the combined fit used in the ^H case to improve 
the fit result. The box in the upper plot shows the same distribu- 
tion with a logarithmic scale to show the ^H component that is not 
visible on a linear scale. The ^H component in the ^H distribution 
in the bottom plot is suppressed by the additional selection shown 
in Fig [6] The hydrogen sample also contains a small fraction of 
•^H events which are used to test the reliability of the Monte Carlo 
simulation. 

way as described in section 13.1.31 for the charge 
misidentification study. The efficiency of the track- 
ing system was however obtained from the Monte 
Carlo simulation. The decrease in efficiency from 
2006 to 2007 is due to the failure of some of the 
front-end chips in the tracking system. This situ- 
ation was incl uded in the Monte Ca rlo simulation, 
as discussed in lAdriani et al. I (J2013D . The efficien- 
cies of the various selections are reported in Table 
1 

• Hadronic interactions: Helium and hydrogen nu- 
clei may be lost due to hadronic interactions in the 
2 mm thick aluminium pressurized container and 
the top scintillator detectors. The correction to the 
flux due to this effect was included in the PAMELA 
geometrical factor as follows: 



G{E) = [1 - b{E)] Gf 



(2) 



where G{E) is the effective geometrical factor used 
for the ffux determination, h{E) is a correction fac- 
tor which accounts for the effect of inelastic scat- 
tering, Gp is the nominal geometrical factor which 
is almost constant above 1 GV and slowly de- 
creases by ~ 2% to lower energies, where the par- 
ticle trajectory in the magnetic field is no longer 



Table 1 

Number of selected and reconstructed events for the Z 



1 sample. 



Kinetic energy 
per nucleon (GeV/n) 



H events 



H events 



H contamination 



^H extrapolated events 
at top of payload (Af-poP ) 



^H extrapolated events 
at top of payload (A''top) 



0.120 - 


-0.132 


38927 ± 477 


2179 ± 99 


207 ±9 


70273 ± 867 


3474 ± 187 


0.132 - 


-0.144 


47689 ± 521 


2401 ± 60 


212 ±5 


85115 ± 940 


4120 ± 106 


0.144 - 


-0.158 


58104 ± 587 


2815 ±114 


230 ±9 


102825 ± 1049 


4800 ± 209 


0.158 - 


-0.173 


70390 ± 635 


3284 ± 120 


248 ±9 


122634 ±1115 


5594 ± 225 


0.173- 


-0.190 


84622 ± 703 


3752 ± 124 


261 ±9 


145937 ± 1227 


6397 ± 227 


0.190 - 


- 0.208 


102408 ± 781 


4312 ± 134 


276 ±9 


177148 ± 1369 


7385 ± 246 


0.208 - 


- 0.228 


123869 ± 843 


4865 ± 144 


286 ±8 


212463 ± 1465 


8343 ± 265 


0.228 - 


- 0.250 


145923 ± 929 


5505 ± 153 


298 ±8 


248196 ± 1601 


9298 ± 278 


0.250 - 


- 0.274 


173016 ± 996 


6015 ± 159 


302 ±8 


291879 ± 1695 


10165 ± 287 


0.274 - 


- 0.300 


200767 ± 1105 


6683 ± 173 


312 ±8 


336058 ± 1868 


11269 ±310 


0.300 - 


- 0.329 


231937 ±1129 


7419 ± 192 


324 ±8 


388972 ± 1903 


12635 ± 345 


0.329 - 


-0.361 


266950 ± 1259 


8274 ± 209 


339 ±9 


443768 ± 2107 


14143 ± 374 


0.361 - 


- 0.395 


303941 ± 1354 


9356 ± 227 


358 ±9 


506860 ± 2291 


15447 ± 402 


0.395 - 


- 0.433 


343790 ± 1389 


10153 ± 245 


358 ±9 


568056 ± 2315 


17005 ± 432 


0.433 - 


- 0.475 


382305 ± 1511 


10922 ± 259 


345 ±8 


630535 ± 2517 


18588 ± 461 


0.475 - 


- 0.520 


424071 ± 1583 


11551 ± 287 


310 ±8 


703994 ± 2654 


20277 ± 511 


0.520 - 


- 0.570 


466939 ± 1658 


12425 ± 337 


258 ±7 


775857 ± 2784 


22210 ± 608 



Table 2 

Number of selected and reconstructed events for the Z 



■ 2 sample. 



Kinetic energy 
per nucleon (GeV/n) 



^He events 



^He events 



*He extrapolated events 
at top of payload (A''top) 



•^He extrapolated events 
at top of payload (A^Top) 



0.126- 


- 0.141 


21092 ± 636 


1522 ± 107 


32094 ± 964 


2387 ± 170 


0.141 - 


- 0.158 


24709 ± 652 


1929 ± 107 


37995 ± 1015 


3023 ± 170 


0.158- 


- 0.177 


29506 ± 695 


2403 ± 131 


45235 ± 1068 


3768 ± 208 


0.177- 


-0.198 


35215 ±817 


2952 ± 143 


54237 ± 1243 


4645 ± 226 


0.198 - 


- 0.222 


41078 ± 897 


3583 ± 162 


64104 ± 1405 


5618 ± 256 


0.222 - 


- 0.249 


47654 ± 980 


4306 ± 163 


74926 ± 1545 


6864 ± 264 


0.249 - 


- 0.279 


53935 ± 1077 


5124 ± 183 


84403 ± 1640 


8203 ± 295 


0.279 - 


- 0.312 


61541 ± 1148 


6040 ± 196 


98143 ± 1851 


9853 ± 325 


0.312- 


- 0.350 


68608 ± 1145 


7046 ± 207 


109058 ± 1837 


11645 ±345 


0.350 - 


- 0.392 


76090 ± 1224 


8125 ± 226 


121874 ± 1971 


13737 ± 388 


0.392 - 


- 0.439 


83651 ± 1268 


9246 ± 255 


135138 ± 2058 


16080 ± 453 


0.439 - 


- 0.492 


90552 ± 1262 


10353 ± 266 


143775 ± 2017 


18371 ± 471 


0.492 - 


-0.551 


99140 ± 1422 


11366 ± 286 


159728 ± 2310 


20577 ± 521 


0.551 - 


- 0.618 


104636 ± 1483 


12173 ±321 


167348 ± 2397 


23111 ±612 


0.618- 


- 0.692 


110854 ±1665 


13209 ± 332 


181468 ± 2771 


26133 ± 662 


0.692 - 


- 0.776 


114097 ±1624 


13946 ± 363 


186092 ± 2665 


28768 ± 762 


0.776 - 


- 0.870 


117080 ± 1768 


13985 ± 362 


193787 ± 2935 


30778 ±811 



Table 3 

Selection efHciencios, divided by detector. 



iH 



^H 



3 He 



^He 



ToF and AC selections 
Tracker selections (2006) 
Tracker selections (2007) 
Tracker dE/dx selection 



88.99% ±0.13% 

88.6% ± 0.2% 

62.2% ± 0.3% 

100% 



88.51% ±0.18% 
87.1% ±0.3% 
60.8% ± 0.4% 

96.67% ±0.13% 



89.0% ± 0.4% 
89.5% ± 0.4% 
61.2% ± 0.6% 
95.5% ± 0.2% 



88.1% ±0.4% 

88.3% ± 0.4% 

62.1% ±0.7% 

100% 



Total efficiency (2006) 
Total efficiency (2007) 



78.8% ± 0.3% 
55.4% ± 0.5% 



74.5% ± 0.5% 
52.0% ± 0.5% 



76.1% ±0.8% 
52.0% ± 0.9% 



77.8% ±0.7% 
54.7% ± 1.2% 



straight. The requirement on the fiducial volume 
corresponds to a geometrical factor Gp ~ 19.9 cm^ 
sr above 1 GV. The correction factor b{E) is dif- 
ferent for each isotope and has been derived from 
the Monte Carlo simulation, being ~ 6% for pro- 
tons, ~ 10% for deuterium, and ~ 13% for both 
helium isotopes. The nominal geometrical factor 
and the effective geometrical factor for each isotope 



are shown in Fig. [7l 

Contamination: The contribution to ^H from ^He 
inelastic scattering was evaluated from the simu- 
lation (Fig. [8]) and subtracted from the raw ^H 
counts (see column 4 of Table [T]). The contami- 
nation in the ^He sample from ^He fragmentation 
was also evaluated and it was estimated to be less 
than 1%. This was included in the systematic un- 




Figure 5. Examples of 1//3 distributions for fielium in the 0.312 - 
0.350 GeV/n kinetic energy range for ^He (top) and ^He (bottom). 
The solid line shows the estimated *He and ^He signal while the 
dashed line shows the combined fit used in the "^He case to improve 
the fit result. 
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Figure 6. The lowest of the 12 energy releases in the tracking sys- 
tem as a function of rigidity for Z = 1 events. Events with energy 
releases above the black line were selected for the ^H analysis. 

certainty of the measurement. 

• Energy loss and resolution: The finite resolution of 
the magnetic spectrometer and particle slowdown 
due to ionization energy losses results in a distor- 
tion of the particle spectr a. A Bayesian un fold- 
ing procedure, described in lD'Agostinil ()1995[ ). was 
used to derive t he number of events a t the top of 
the payload fsee lAdriani et aU\ ()2011j )). 

The flux was then calculated as follows: 

*^°P(^) " TG{E)AE ^^^ 

where TVtoP (E) is the unfolded particle count for energy 



E, also corrected for all the selection efficiencies (see Ta- 
bles [1] and [2l rightmost two columns), AE is the energy 
bin width, and G{E) is the effective geometrical factor. 
The live time, T, as evaluated by the trigger system, 
depends on t he orb ital selection as described in section 
I3.1.2l fe.g. see [Bruno (2008)). The live and dead time are 
cross-checked with the total acquisition time measured 
by the on-board CPU to remove possible systematic ef- 
fects in the time counting. 

3.4. Systematic uncertainties 

The possible sources of systematic uncertainties con- 
sidered in this analysis arc listed below and are also in- 
cluded in Tables Hand El and in Figs. H [TO] and [111 

• Quality of the 1//3 fit: The quality of the Gaus- 
sian fit procedure was tested using the truncated 
mean of the energy deposited in the electromag- 
netic calorimeter to select pure samples of ^H and 
^H from non-interacting events. The two sam- 
ples were then merged to form a control sample 
for the fitting algorithm. The number of recon- 
structed events from the Gaussian fit was found 
to agree with the number of events selected with 
the calorimeter, so no systematic uncertainty was 
assigned to this procedure. 

• Selection efficiencies: The estimation of the selec- 
tion efficiencies is affected by a statistical error due 
to the finite size of the sample used for the effi- 
ciency evaluation. This error was considered and 
propagated as a systematic uncertainty. For the 
efficiency of the ToF and AC selections this un- 
certainty is 0.21% at low energy (120 MeV/n) and 
drops to 0.14% at high energy (600-900 MeV/n). 
For the tracker selections the uncertainty is 0.3% 
at low energy increasing to 0.4% at high energy. 

• Galactic particle selection: The correction for par- 
ticles lost due to this selection has an uncertainty, 
due to the size of the Monte Carlo sample, which 
decreases from 6% to 0.07% as the energy increases 
from 120 MeV/n to 900 MeV/n. 

• Contamination subtraction: The subtraction of the 
contamination results in a systematic uncertainty 
on the ^H flux of 1.9% at low energy dropping be- 
low 0.1% at 300 MeV/n due to the finite size of the 
Monte Carlo sample. To test the validity of the 
Monte Carlo simulation the '^H component, identi- 
fied as the additional cluster of events at low f3 in 
the hydrogen sample visible in Fig. [4l was used. 
The •^H events arc created by '^Hc spallation in the 
top part of the apparatus since no tritium of galac- 
tic origin should survive propagation to Earth. The 
observed number of ^H events was used to test that 
the Monte Carlo simulation correctly inferred the 
number of ^H events coming from *He fragmen- 
tation. For example, for the 2006 data-set in the 
rigidity range between 1.7 GV and 1.8 GV the flight 
data sample contains 136 ± 17 tritium events, while 
110±15 events are expected according to the Monte 
Carlo simulation (Fig. [4]). Simulation and flight 
data were in agreement within a 10% tolerance. 




Figure 7. The nominal geometrical factor Gp as a function of 
rigidity (solid line). The filled bands represent the effective geo- 
metrical factor G{E) for each isotope with the associated uncer- 
tainty. 
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Figure 8. Raw counts for ^H collected by PAMELA in 2006 to- 
gether with the expected contamination from He spallation in the 
top part of the apparatus. 

This discrepancy was treated as an additional sys- 
tematic uncertainty on the estimated number of 
contamination events. The 10% systematic uncer- 
tainty on the ^H contamination translates in an 
additional 1% uncertainty on the number of recon- 
structed ^H events. 

• Geometrical factor: The uncertainty on the effec- 
tive geometrical factor as estimated from the Monte 
Carlo simulation is almost independent of energy 
and amounts to 0.18%. 

• Unfoldin g p ro cedure : As discussed in 
lAdriani et al. I (|2011f ) two possible systematic 
effects have been studied regarding the unfolding 
procedure: the uncertainty associated with the 
simulated smearing matrix and the intrinsic 
accuracy of the procedure. The former was con- 
strained by the checking for compatibility between 
measured and simulated spatial residuals and was 
found to be negligible. The latter was estimated 
by folding and unfolding a known spectral shape 
with the spectrometer response and was found to 
be 2%, independent of energy. 



4. RESULTS 

Figure [9] and [10] show hydrogen and helium iso- 
tope fluxes (top) and the ratios of the fluxes (bot- 
tom). Results are also reported in Tables [J] and 
[5l Fig. [TT] shows the ^H/'*He ratio as a function 



of kinetic ener gy per nucleon, compared to previous 
measurements (lAguilar et al. 1 120111: iWang et al. 1 120021: 
Reimcr et aT\ 119981: IWefel et al. I 119951: iWebber et al. I 
1991: Bcat tv et al. 1119931 ). It is worth noting that flux ra- 
tios (in particular the "^He/^He ratio), if modulated usin g 
the force-field approximation (jGleeson fc Axfordlll968D . 
show very little dependence on solar activity and can 
therefore be used to discriminate between various prop- 
agation models of GCR in the Galaxy. 

The PAMELA results are the most precise to date. Con- 
sidering the relatively large spread in the existing data, 
PAMELA results agree with previous measureinents, in 
particular with BESS results for ^H and IMAX results 
for '^He. Previous measurements are affected by large 
uncertainties and, for '^He where more measurements are 
available, there is a large spread between data. All the 
measurements displayed in Figures [9l I10[ and lll[ except 
AMS-01, are from balloon-borne experiments and are af- 
fected by a non-negligible background of atmospheric sec- 
ondary particles. 

A high precision measurement of the H and He iso- 
tope quartet abundances represents a significant step for- 
ward in modelling the origin and propagation of GCRs. 
The constraints on diffusion-model parameters set by 
the qua rtet (^H, ^H, ^He , and ^He) were recently re- 
visited (jCoste et al. 11201^ . It was found that the con- 
straints on the parameters were competitive with those 
obtained from the B/C fiux ratio analysis and available 
data supported the universality of GCR propagation in 
the Galaxy. The tightest constraint was obtained when 
the He flux was included in the fit. This is because at 
energies of a few GeV about 10% of He is from fragmenta- 
tion of heavier nuclei, which is a non-negligible amount 
given the precision (1% s tatistical) of the PAMELA He 
data (|Adriani e^ al. 11201 1[ ). 
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Table 4 

Hydrogen isotope fluxes and their ratio, errors are statistical and systematics respectively. 



Kinetic energy 

at top of payload 

(GeV n-l) 



iR flux 



(GeV n- 



s sr) 



^H flux 



(GeV n ^ m-^ s sr) 



^H/ iR 



0.120 


0.132 


(9.86 ±0.15 ±0.31) lO^ 


(30.8 ±1.8 ±1.0) 


(3.12 ±0.19 ±0.20) 


10- 


-2 


0.132 


0.144 


1.064 ±0.014 ±0.030) 


10=* 


(33.0 ±1.6 ±0.9) 


(3.10 ±0.15 ±0.18) 


10- 


_2 


0.144 


0.158 


1.136 ±0.013 ±0.029) 


10=* 


(34.2 ±1.4 ±0.9) 


(3.01 ±0.13 ±0.16) 


10- 


_2 


0.158 


0.173 


1.182 ±0.012 ±0.028) 


10=* 


(35.1 ±1.4 ±0.8) 


(2.97 ±0.12 ±0.14) 


10- 


_2 


0.173 


0.190 


1.233 ±0.012 ±0.027) 


10=* 


(35.3 ±1.3 ±0.8) 


(2.86 ±0.11 ±0.13) 


10- 


_2 


0.190 


0.208 


1.300 ±0.011 ±0.027) 


10=* 


(36.6 ±1.3 ±0.8) 


(2.81 ±0.10 ±0.12) 


10- 


-2 


0.208 


0.228 


1.351 ±0.010 ±0.027) 


10=* 


(37.1 ±1.2 ±0.7) 


(2.74 ±0.09 ±0.11) 


10- 


-2 


0.228 


0.250 


1.405 ±0.010 ±0.027) 


10=* 


(36.3 ±1.1 ±0.7) 


(2.58 ±0.08 ±0.10) 


10- 


-2 


0.250 


0.274 


1.453 ±0.009 ±0.027) 


10=* 


(34.5 ±1.0 ±0.6) 


(2.37 ±0.07 ±0.09) 


10- 


-2 


0.274 


0.300 


1.486 ±0.008 ±0.027) 


10^ 


(33.9 ±0.9 ±0.6) 


(2.28 ± 0.06 ± 0.08) 


10- 


-2 


0.300 


0.329 


1.514 ±0.008 ±0.027) 


10=* 


(34.1 ±0.9 ±0.6) 


(2.25 ±0.06 ±0.08) 


10- 


_2 


0.329 


0.361 


1.517 ±0.007 ±0.026) 


10=* 


(33.6 ±0.9 ±0.6) 


(2.22 ±0.06 ±0.08) 


10- 


_2 


0.361 


0.395 


1.520 ±0.007 ±0.026) 


10=* 


(32.0 ±0.8 ±0.6) 


(2.10 ±0.06 ±0.07) 


10- 


-2 


0.395 


0.433 


1.484 ±0.006 ±0.025) 


10=* 


(31.4 ±0.8 ±0.6) 


(2.11 ±0.05 ±0.07) 


10- 


_2 


0.433 


0.475 


1.480 ±0.006 ±0.025) 


10=* 


(30.6 ±0.8 ±0.5) 


(2.07 ±0.05 ±0.07) 


10- 


_2 


0.475 


0.520 


1.499 ±0.006 ±0.025) 


10=* 


(29.9 ±0.8 ±0.5) 


(2.00 ±0.05 ±0.07) 


10- 


_2 


0.520 


0.570 


1.429 ±0.005 ±0.024) 


10=* 


(29.0 ±0.8 ±0.5) 


(2.03 ±0.06 ±0.07) 


10- 


-2 



Table 5 

Helium isotope fluxes and their ratio, errors are statistical and systematics respectively. 



Kinetic 


energy 


^He flux 




=*Hc flux 


=5He/ 4Hc 




at top of 


payload 












(GeV 


n-l) 


(GeV n-l m2 s sr) 


-1 


(GeV n-l m2 s sr)"! 






0.126 - 


0.141 


2.03 ± 0.06 ± 0.04) 


102 


(17.1 ±1.8 ±0.4) 


(8.4 ±0.9 ±0.4) -10-2 


0.141 - 


0.158 


2.20 ±0.06 ±0.04) 


10^ 


(21.3 ±1.7 ±0.5) 


(9.7 ±0.8 ±0.4) -10-2 


0.158 - 


0.177 


2.27 ±0.06 ±0.04) 


10^ 


(23.4 ±1.5 ±0.5) 


(1.03 ±0.07 ±0.04) 


10-i 


0.177 - 


0.198 


2.33 ±0.06 ±0.04) 


10^ 


(24.8 ±1.4 ±0.5) 


(1.06 ±0.06 ±0.04) 


10-i 


0.198 - 


0.222 


2.38 ±0.05 ±0.04) 


10^ 


(25.4 ±1.2 ±0.5) 


(1.07 ±0.06 ±0.04) 


10-1 


0.222 - 


0.249 


2.42 ±0.05 ±0.04) 


10^ 


(26.7 ±1.2 ±0.5) 


(1.10 ±0.05 ±0.04) 


10-1 


0.249 - 


0.279 


2.39 ±0.05 ±0.04) 


10^ 


(27.3 ±1.0 ±0.5) 


(1.14 ±0.05 ±0.04) 


10-1 


0.279 - 


0.312 


2.37 ±0.04 ±0.04) 


10^ 


(28.2 ±1.0 ±0.5) 


(1.19 ±0.05 ±0.04) 


10-1 


0.312 - 


0.350 


2.29 ±0.04 ±0.04) 


10^ 


(28.4 ±0.9 ±0.5) 


(1.24 ±0.04 ±0.04) 


10-1 


0.350 - 


0.392 


2.22 ±0.04 ±0.04) 


10^ 


(28.9 ±0.8 ±0.6) 


(1.30 ±0.04 ±0.05) 


10-1 


0.392 - 


0.439 


2.11 ±0.03 ±0.04) 


10^ 


(29.2 ±0.8 ±0.6) 


(1.39 ±0.04 ±0.05) 


10-1 


0.439 - 


0.492 


1.93 ±0.03 ±0.03) 


10^ 


(28.6 ±0.8 ±0.6) 


(1.49 ±0.04 ±0.06) 


10-1 


0.492 - 


0.551 


1.84 ±0.03 ±0.03) 


10^ 


(27.3 ±0.7 ±0.6) 


(1.49 ±0.04 ±0.06) 


10-1 


0.551 - 


0.618 


1.72 ±0.02 ±0.03) 


10^ 


(26.5 ±0.7 ±0.5) 


(1.55 ±0.05 ±0.06) 


10-1 


0.618 - 


0.692 


1.59 ±0.02 ±0.03) 


10^ 


(25.4 ±0.7 ±0.6) 


(1.60 ±0.05 ±0.06) 


10-1 


0.692 - 


0.776 


1.44 ±0.02 ±0.03) 


10^ 


(24.5 ±0.6 ±0.5) 


(1.70 ±0.05 ±0.07) 


10-1 


0.776 - 


0.870 (1 


292 ± 0.020 ± 0.023 


)-102 


(23.0 ±0.6 ±0.5) 


(1.78 ±0.05 ±0.07) 


10-1 
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Figure 9. ^H and ^H a bsolute fluxes (a) and their ratio (b) compared to previous experiments: AMS-01 I IAguilar et al. ||2002I . [20111 ). 
BESS (Wa nir~ei al. II2002I ). IMAX IjReimer et al. iri998) . Error bars show the statistical uncertainty while shaded areas show the systematic 
uncertainty. 
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Figure 10. ^Hc and '^Hc absolute fluxes (a) and their ratio ( b) compared to previ ous experiments: AM S l lAguilar et al. | |2011l). BES S 
HWang et al. 2002), IMAX (Reimer et ai. I9m\. SMILI-2(Wef cl et al. |[T995l 'l. MASS (Webber et al. II19911 '|. SMILI-1 IjBeattv et al II1993I 'I. 
Error bars show statistical uncertainty while shaded areas show systematic uncertainty. 
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Figure 11. 'H/'*H e ratio compared to previous experiments: AMS-01 l lAguilar et al. I l20Tll) . BESS llWang et al. I l200l) . IMAX 
HReimer et al. iri998l ). Error bars show statistical uncertainty while shaded areas show systematic uncertainty. 



